% 创建时间：2025-06-20 23:11:31
% 作者：兪馳
clc;
clear;

A=[1 2 3;
    4 5 6;
    7 8 0]; % 定义矩阵 A
C=-[10,5,4;
    5,6,7;
    4,7,9]; % 定义矩阵 C
X=lyap(A,C), % 使用 lyap 函数求解 Lyapunov 方程 AX + XA' + C = 0 的解 X
norm(A*X + X*A' + C) % 计算 AX + XA' + C 的范数，检验解的准确性